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Abstract 

Chaotic systems can be synchronized by linking them to a common signal, subject to certain condi¬ 
tions. However, the presence of multiple driving signals coming from different systems, give rise to novel 
behavior. The particular case of Lorenz systems, with two independent systems driving another system 
through drive-response coupling has been studied in this paper. This is the simplest arrangement which 
shows the effect of “frustrated synchronization” due to competition between the two driver systems. The 
resulting response system attractor deviates significantly from the conventional Lorenz attractor. A new 
measure of desynchronization is proposed, which shows a power-law scaling relation with the competition 
parameter. 
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1 Introduction 


The synchronization of chaotic systems is a difficult problem owing to their extremely sensi¬ 
tive dependence on initial conditions. Any initial correlation present between identical systems, 
starting from very close initial conditions, exponentially decrease to zero with time. Thus, for 
all practical purposes, any initial synchronization between the systems is bound to disappear 
rapidly. In recent times, however, some methods of achieving synchronized behavior between 
chaotic systems have been proposed. Pioneering work in this respect has been done by Pecora 
and Carroll |]I|, who used the concept of a response system locking on to a driver system. So far, 
such studies have been limited to driving a response system by a single driver system. However, 
the knowledge gained from studying such simple systems may not be adequate to give us an idea 
as to how systems consisting of multiple independent driver systems, competing with each other 
to synchronize the same response system, will behave. The Pecora-Carroll driving mechanism 
can be seen as the “strong-coupling” limit of a general scheme of directionally- oriented couplings 
in a network of chaotic elements. 


The synchronization of bidirectionally coupled chaotic systems is stable provided the coupling 
strength is at least half the Lyapunov exponent of the system Q . One-way coupling (or, driving 
one chaotic system by another) can also lead to synchronization, provided certain conditions 
are satisfied |]|], |^, [Q]. The drive-response method consist of the following steps. First an 
n-dimensional autonomous system 


dx 

dt 


F(x), 


is divided into two parts, driving (x^) and responding (x,,): 


dxrf 

dt 


g(xd,x.r). 


dxr 

dt 


h(xrf,Xr), 


where, Xrf = {xi,.. .,Xm), g = [/i(x),...,/^(x)], x,, = (x^+i,..., and h = [/m+i(x),...,/„(x)]. 
A replica subsystem x). identical to x^ is then created and driven with the x^ variables of the 
original system. Therefore, the replica subsystem equations are. 


dx), 

dt 


h(xrf, x(.). 


The responding subsystems x^ and x). will synchronize only if dx,, = |xj, — x(,| —0. According 
to Pecora and Carroll, this occurs if and only if the conditional Lyapunov exponents of the x,. 
subsystem are all negative. 


Drive-response synchronization has been realized in various electrical circuit experiments. It has 
also been used in experiments of secure communication where a chaotic masking signal is added 
to the transmitted signal. It is then recovered at the receiving end by subtracting the chaotic 
signal regenerated by synchronization Q. 
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Besides the Pecora-Carroll method, other synchronization procedures have also been proposed. 
Of these, the Variable Control Feedback (VCF) method is of particular interest, as it can be 
used for both control and synchronization of chaos |^]. In fact, the Pecora-Carroll method turns 
out to be a special limiting case of this method. VCF consists of adding a feedback term to 
a dynamical system to guide it into some prescribed state. If ^ = F(x) be an n-dimensional 
dynamical system and x* be the desired state to which the system has to be brought, then VCF 
involves modifying the system dynamics to: 

|=F(x)-A(x-x-) 

where A is the set of n feedback multipliers. If x* be the output of a chaotic system F'(x), then the 
system synchronizes with F(x). In the large-A limit, VCF reduces to the Pecora-Carroll method. 
Specihcally, the feedback parameters for the driving subsystem variables, A^ —> oo, while the 
remaining As are set to zero. 

In this paper some observations have been reported on the attractor structure of a chaotic system 
which has been subjected to simultaneous synchronization by two other identical chaotic systems 
competing with each other. Section 2 introduces the model used for studying competition among 
synchronizing chaotic systems and includes a short analysis of the fixed points and their stability. 
Section 3 contains the results of computer simulations of the system. Finally, possible directions 
of future research and the relevance of this type of research to the theory of neural computation 
are discussed. 


2 Competition among synchronizing Lorenz systems 


The investigation of competition among synchronizing chaotic systems was carried out using the 
Lorenz system of equations [^ , Q . This well-known paradigm of chaos is defined by the following 
set of equations: 


dx 

dt 

dt 

dz 

dt 


= a {y - x), 
rx — y — xz, 


= xy — bz, 


( 1 ) 

( 2 ) 

( 3 ) 


where, a, r and b are real, positive parameters. There are three fixed points for this system: 
Fi = (0,0,0), F2 = - 1), 'JK'f - 1)> ^ - 1); and F3 = ( -sjb{r - 1), -^/b{r^^, r - 1). 

The local stability of the fixed point {xj , yj , Zf) is determined by the eigenvalues of the 
Jacobian 

— a a 0 


{ r — Zf ) — 1 — Xf 

Vf Xf -b 


( 4 ) 
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Evaluation of the matrix shows that for 0 <r<l,Fiis the only stable fixed point. For r > 1 , Fi 
becomes unstable and the phase-space trajectory of the system converges to either F2 or F3. For 
r > Vc = a{a + b + 3 )/{a — b — 1 ) the system’s trajectory perpetually wanders along the extremely 
complicated structure of the stable and unstable manifolds of the fixed points, exhibiting chaotic 
behavior. 


For the present work the effect of two driving systems, designated as driving systems 1 (xi, yi, zi) 
and 2 {x2,y2, Z2), competing to synchronize a responding system {x3,y3,Z3) was studied. The 
responding system was driven using the y variable. A competition parameter a was defined to 
indicate the strength of the driving systems relative to each other. The maximum value of a was 
normalized to unity. Therefore, the y variable of the responding system was defined in terms of 
the two driving systems as: 

2/3 = ayi + {l- a)y 2 . ( 5 ) 

We consider first the case where the two driving systems have the same r-parameter value, and 
then, the more general case, where the two r-values are different (ri and r2, say). The a and 
6-parameter values are considered to be the same in all cases. 

Case I: ri = r2 = r 


It is obvious that for a = 1 the responding system synchronizes with driver system 1 , whereas for 
a = 0 , it synchronizes with system 2 . The attractor of the response system, is identical to that 
of the conventional Lorenz system (fig. 1 (a)). For 0 < a < 1 , the responding system (^3,^3, 2:3) 
has nine fixed points: 

Fi = (0,0,0), 

F2 = (A/6(r-l),A/6(r-l),r-l), 

F 3 = 1), 

F4 = {a^/b{r - l),a^/b{r - 1), a'^{r - 1)), 

F5 = ((1 - a)^yb{r - 1), (1 - a)y^b{r - 1), (1 - af{r - 1)), 

Fe = {-ay^b(r^^,-ay^bir^^,a‘^{r - 1)), 

F7 = (-(l-a)A/ 6 (f^Fiy,_(i_a)y 5 (fWi), (i_a) 2 (^_l)), 

Fs = ((2a - l)\/6(r - 1), (2a - l)y'b{r - 1), (2a - lf{r - 1)), 

Fg = (-( 2 a - 1 )V6(r - 1 ), -( 2 a - l)^b{r - 1 ), ( 2 a - - 1 )). 


Note that the first three fixed points are those of the uncoupled Lorenz system. To find out 
about the stability of these fixed points we need to calculate the eigenvalues of the corresponding 
Jacobian, Jh The partially block-diagonal form of the matrix makes the calculation easy: 


J' 


J 03x3 03x2 

03x3 J 03x2 ; 

A B 


( 6 ) 


where, J is the Jacobian (eqn. 4 ) of the unperturbed Lorenz system of equations, O^xn is a null 
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matrix having m rows and n columns, and the other matrices are defined as, 


, 0 a a 0 

A = 

0 a Xf^ 0 

g ^ 0 (1 - a)a 0 

0 ( 1 - a ) x/3 0 

and, 

3r = ~ ° 

a y/i + ( 1 - a ) 2//2 - b 

Here fk refers to the fixed point of the kth Lorenz system. 


(7) 

( 8 ) 

( 9 ) 


For 0 < r < 1 , the only stable fixed point is Fi. For r > 1 , Fi loses its stability, and there are 
four new stable fixed points: F2,F3,F8 and Fg. For r > Vc = (j{a + 6 + 3 )/{a — 6 — 1 ), these 
fixed points lose their stability and the system shows only chaotic behavior. The most interesting 
instance is that of a = 0 . 5 , where maximal competition occurs. In this case, Fg = Fg = Fi, 
F4 = F5 and Fg = F7 (fig. 2 ). The attractor of the responding system is found to be stretched 
over its 3 -dimensional phase space showing an extremely tangled structure (fig. 1 (b)). This is 
due to the extremely complicated motion of the response system trajectory along the stable and 
unstable manifolds of the fixed points Fi,F2,F3,F4 and Fg. The coupling with driver system 1 
tries to force the response system into synchronization with it, but at the same time, the coupling 
with driver system 2 desynchronizes the trajectory. The synchronization is therefore ‘frustrated’ 
by the competition between the two driver systems. The “frustrated” response system attractor 
reduces to the conventional Lorenz attractor if a —> 0 or 1 , when competition is absent. 

The attractor structure is found to be quite robust. If we start from two different initial conditions 
for the responding system, {x,y,z) and {x',y',z'), say, then for stable synchronization, the two 
respective trajectories should converge rapidly. However, whereas in the Pecora-Carroll case, 
convergence occurs to the standard Lorenz attractor, in this case, both the trajectories converge 
to the “frustrated” attractor. 


The stability of synchronization can be demonstrated analytically by linear stability analysis of 
the error dynamics. Defining the dynamical error between two response system trajectories (x 
and x') which have different initial conditions, as e = x — x^, the error equations can be written 
as: 


dCr 

(10) 

€y 0, 

(11) 

{ayi -h (1 - a)y2)ex - be^- 

(12) 


Here we have assumed that the equation parameters for the two systems are identical. The 
error system of equations has an equilibrium point at e = (0, 0, 0), which corresponds to perfect 


5 



synchronization. The local stability of synchronization can then be checked by looking at the 
eigenvalues of the Jacobian of the error equations; 




— a 0 

ayi + (l-a)y2 -h 


( 13 ) 


The eigenvalues are —a and —6, which are the conditional Lyapunov exponents of the response 
system. As both eigenvalues are negative, the synchronization is locally stable, and any difference 
in initial conditions rapidly goes to zero. Note that, this does not prove the global stability of 
the synchronized state. However, simulations have verified that even in the presence of large 
deviations in initial conditions, synchronization with the “frustrated” trajectory is achieved. 
This indicates that, although exact synchronization with the driver system cannot be achieved, 
the “frustrated” system can still be used for secure communication through chaotic masking. 
This has been established through simulations reported below. 


Case II: ri ^ r2 


When the value of the r-parameter of the two driving systems is not the same, the fixed points 
are given by; 

Fi = ( 0 , 0 , 0 ), 

F2 = {a^Jb{ri - 1) + (1 - a)y'b{r 2 - 1), a^b{r - 1) + (1 - a)y'b{r 2 - 1), - 1) + (1 - af(r 2 - 

1) + 2a(l - a)^/lrl - l)(r2 - 1)), 

F3 = {-a^/b(r^'^^ - (1 - a)yjb{r2 - 1 ), -a^Jb{r - 1 ) - (1 - a)^Jb{r2 - 1 ), a^(ri - 1 ) + (1 - 
a)2(r2 - 1) + 2a(l - a)^{n - l)(r2 - 1)), 

F4 = {a\/b{ri - l),a^Jb{ri - - 1)), 

F5 = ((1 - a)yjb{r 2 - 1), (1 - a)^/b{r 2 - 1), (1 - a)‘^{r 2 - 1)), 

Fg = {-ay/b{ri - 1), -a^b{ri - - 1)), 

Fy = (-(1 - a)^/b{r2 - 1), -(1 - a)^b{r2 - 1), (1 - a)2(r2 - 1)), 

Fs = ia^yb{rl - 1 ) - { 1 -a)^b{r 2 - l),a^yb{r - 1 )-{l-a)^yb{r 2 - 1 ), {n - 1 ) + (1 -a)'^{r 2 - 
1) - 2a(l - a)vTn - 1)(?’2 - 1)), 

Fg = {-ay^b{ri - 1 ) + (1 - a)^b{r2 - 1 ), -a^b{r - 1 ) + (1 - a)^b{r2 - l),a'^{ri - 1 ) + (1 - 
af(r2 - 1) - 2a(l - a)\/(n - l)(r2 - 1)). 

Fig. 3 shows the (ri, r2)-parameter space. The stable fixed points at different regions are indicated 
in the diagram. The dotted line corresponds to the special case ri = r2 which has been considered 
above. Note that, whereas in the general case all the fixed points are stable in some region or 
other, in the special case of ri = r2, four of the fixed points, viz., F4,F5,F6 and Fy, are always 
unstable. When one of the r-values go over to the chaotic regime, while the other r-value remains 
fairly below it, asymptotic synchronization with the chaotic trajectory is observed Q. The time 
required to ultimately synchronize with the chaotic attractor is a function of both the r-parameter 
values. The synchronization is phase- synchronization rather than state- synchronization, as the 
response system chaotic attractor is a scaled replica of the driver system attractor. The scaling 
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factor is a for synchronization with driving system 1, and (1 — a), for driving system 2. When 
both the r-values are in the chaotic regime, the “frustrated synchronization” situation occurs. 


3 Simulation Results 


For conducting simulations, the parameter values chosen were ri = r 2 = 28, a = 10 and b = 8/3. 
The trace of the Jacobian (which is equal to the sum of the Lyapunov exponents) for the total 
system, including the driver and response systems, is -40.0. So the overall system is diffusive 
and possesses an attractor. The competition parameter a was varied in the interval [0,1]. The 
differential equations were numerically solved using the fourth-order Runge-Kutta method with 
step-size = 0.025. The phase-space trajectory of the responding system {x 3 ,y 3 ,Z 3 ) was observed 
with different values of a from t = 0 to t = 100 . At the limit a = 0 (or 1 ) the responding 
system trajectory is identical to that of a unperturbed Lorenz system (fig. 1(a)). However, as 
a 0.5 (where maximal competition occurs), the trajectory deviates more and more from the 
standard Lorenz form. At a = 0.5, the trajectory moves in a complicated path around the fixed 
points Fi, F 2 and F 3 (note that, at a = 0.5, Fg = Fg = Fi (fig. 1(b)). It appears that for a=0.5, 
the z-variable time-series is much more correlated. This becomes clearer on taking a Fourier 
transform of the data. The power spectral density of the frustrated attractor time-series is low 
in the high-frequency end compared to the unperturbed system time-series. 

The Lyapunov exponents were calculated using Gram-Schmidt technique to create an orthonor¬ 
mal basis every 0.5 seconds of simulation time (this time interval being roughly half the “period” 
of the Lorenz system) and then averaging over 100 iterations. As expected, of the eight expo¬ 
nents, six correspond to those for the two unperturbed driving Lorenz systems (0.84, 0, -14.51). 
The remaining two exponents are the conditional Lyapunov exponents of the responding system 
: -8/3 and -10. This implies the robustness of the “frustrated” attractor - as any deviation from 
the attractor rapidly diminishes. 

To study the degree of synchronization, ^-coordinates of the responding system state (z 3 )were 
plotted against the ^-coordinates of each of the driver system states (zi,Z 2 ), for different values 
of a. If the two are synchronized, the plot gives a straight line. This suggests that the linear 
correlation coefficients, r, between the driver and response system time series, can be used to 
obtain a quantitative measure of synchronization. The linear correlation coefficient between two 
time series data x{t) and = 1 ,..., n), is given by 

^ _ 1 Er=i(a;(i)-x)( 2 /(i)-y) 

' x,y — 5 

n axdy 

where x and ax are the mean and standard deviation respectively, for the time series x{t). A 
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measure of desynchronization is defined as 


<5—1 fz2,Z2,- (14) 

At 0=0, where there is exact synchronization between driver system 2 and the response system, 
(5 = 0. This is a particularly robust measure, as (5 —> 0 for both state- and phase- synchronization. 
The variation of 5 with a is shown in a logarithmic plot (fig. 4). The linear nature of the curve 
over at least 3 orders of magnitude as o —> 0, indicates the presence of a power-law scaling relation 
of the form; 

5 ~ o^, (15) 

where the scaling exponent, /3 ~ 2.0. The scaling exponent was also obtained for r= 50 and 70. 
In both cases, (3 ~ 2.0 within simulation error. 


Another interesting feature studied was the fractal correlation dimension of the frustrated at¬ 


tractor (fig. 5), calculated using the FD3 (ver. 0.3) software ||l^. For the unperturbed Lorenz 
system, this is very close to 2, as the attractor is almost 2-dimensional. As a increases from 
0 to 0.5, the attractor deviates from this two-dimensional shape, which can be quantitatively 
measured by the correlation dimension. As a —>■ 0.5, the attractor structure stretches out more 
and more over the three-dimensional space. This type of enhanced diffusion in phase space seems 
to be a generic feature of frustration in chaotic systems, and has been reported previously in the 
case of Coupled Map Lattices |n|. 


The simulations also showed the robustness of the “frustrated” attractor. Starting from different 
initial conditions, the response system trajectory was found to converge to the same attractor 
structure. This indicates that even in the absence of exact synchronization with any of the 
driver systems, the response system trajectory can be used as a chaotic masking signal for secure 
communication |^. This was verified by adding a small amplitude periodic signal (e.g., a sine 
wave of frequency to = 1/200) to the response system y-variable time series. The resultant time 
series appears to be devoid of any periodic component (fig. 6, top). It is then used to drive another 
Lorenz system, and the x-variable time series of the two systems are subtracted from each other 
to retrieve the original signal (fig. 6, bottom). The modulation of the competition parameter,a, 
by a binary signal for chaotic switching, is another possibility of using the competitive scheme 
for secure communication. 


4 Discussion 


The competitive scheme described here for y-variable coupling was also implemented for x- and 
z-coupling of Lorenz systems. In the former, similar generalized attractor structure was observed, 
while in the latter, where the Pecora-Carroll synchronization does not work, no such structure 
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could be observed. The work done here on coupled Lorenz systems can be extended to other 
systems defined by autonomous set of differential equations as well as discrete maps. However, 
it might be interesting to consider the result of competition in synchronizing non-autonomous 
systems (e.g., the Duffing oscillator). As such systems already have a forcing term present, 
which brings about the onset of chaos, the introduction of additional forcing terms can lead to 
qualitatively new behavior. 


Competitive synchronization in extended systems might also lead to interesting phenomena. Lat¬ 
tices of (globally or diffusively) coupled chaotic elements, where each element can be used both to 
drive other elements, as well as respond to driving signals from yet another set of elements, and 
hence by a series of feedbacks drive its own driving systems, will serve to illustrate interactions be¬ 
tween multiple competing synchronizing feedback loops. The motivation for such a study is that, 
in the human brain, synchronization of activity among different neurons appear to have an impor¬ 
tant functional role in the proper performance of perceptual tasks. It is to be noted that, single 
neurons are capable of chaotic behavior. As the brain is composed of densely connected networks 
of neurons, there is bound to be competitive synchronizing interactions between neural assemblies 
[12|. A dynamic competition parameter, which causes synchronization-desynchronization transi¬ 
tions between various neural sub-assemblies, is a possible mechanism for information processing 
in biological systems. The resultant dynamics will be radically different from the one we are led 
to expect by observing the dynamics of single neurons or small groups of neurons. 


The above work describes the simplest competitive scenario which can show a qualitatively dif¬ 
ferent dynamics from that in the non-competitive situation. It is at present not known how the 
nature of synchronization and the attractor structure of the responding system might be altered 
by increasing the number of competing driver systems. In the brain, where each neuron is con¬ 
nected to ~ 10“^ other neurons, the competitive situation is bound to be far more complicated. 
The manner in which such an extremely competitive synchronization scenario might influence the 
way in which neural networks perform computations and process information is a very interesting 
problem for the future. 
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Fig. 1 The response system attractor for (a) a = 1.0 and (b) a = 0.5 (ri = r 2 = 28, cr = 10, b = 
8/3). 

Fig. 2 The z-coordinate of fixed points of the response system for 0 < a < 1. 

Fig. 3 The (ri, r 2 )-parameter space showing the stable fixed points of the response system at 
different regions. 

Fig. 4 Log-scale plot of desynchronization (<5) for 0 < a < 1. The power-law scaling relation 
(with characteristic exponent, (3 ~ 2.0) is indicated by the solid line fitted to the simulation data. 

Fig. 5 Correlation dimension of the response system attractor for 0 < a < 1. 

Fig. 6 Chaotic masking: the x-variable time series of response system (top); the periodic signal 
obtained by subtracting the regenerated time series from the chaotic carrier wave (bottom). 
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